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Abstract —In this paper, we study the Simultaneous Wireless 
Information and Power Transfer (SWIPT) in a Single-Input 
Single-Output (SISO) two-user Orthogonal Frequency Division 
Multiplexing (OFDM) Interference Channel (IFC). We assume 
that the transmitters are non-cooperative and have perfect 
knowledge of the local Channel State Information (CSI). We 
show that the necessary condition for the optimal transmission 
strategy at high SNR is for the energy transmitter to transmit 
its signal hy allocating its transmit power on a single suhcarrier. 
Accordingly, we propose a one-suhcarrler selection method for 
the energy transmitter and identify the achievable rate-energy 
region. In addition, we further enlarge the achievable rate-energy 
region by enabling a basic form of transmitter cooperation where 
messages are exchanged to inform the energy transmitter about 
the subcarriers unutilized by the Information transmitter. 

I. Introduction 

Recently, Radio Frequency (RE) has been considered as a 
potential resource for energy harvesting in wireless system. 
Because an RE signal carries both information and power, 
SWIPT has been investigated [l]-[8]. There have been several 
studies on SWIPT addressing interference channel (IFC) ||T)- 
0. SWIPT in a two-user and K-user MIMO system have 
been considered in |[T| and 0, respectively. 0 has considered 
SWIPT with partial CSI at the transmitter in MIMO IFC. 
SWIPT in a multi-user MISO IFC has been discussed in lH. 
0 has considered SWIPT in a two-user and K-user SISO IFC. 
Papers 0, 0 have discussed SWIPT in OFDM system but 
did not address the interference channel. 

In this paper, we leverage our past results on SWIPT in 
MIMO IFC 111-0 to address SWIPT in SISO-OFDM IFC. 
To the best of the authors’ knowledge, SWIPT for an OFDM 
interference channel has not been addressed yef We show 
that the optimal non-cooperative transmission strategy at the 
energy transmitter with local CSI knowledge at the transmitter 
(CSI of the links between a transmitter and all receivers) at 
high Signal-to-Noise Ratio (SNR) is to allocate its transmit 
power to a single subcarrier. In this paper, the optimality is 
at high SNR under the assumption that the transmitters are 
non-cooperative and have perfect knowledge of their local 
CSI. Accordingly, we propose a single subcarrier transmission 
startegy and identify the achievable rate-energy region. In 
addition, we show that by enabling some form of transmitter 
cooperation where messages are exchanged to inform the 
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energy transmitter about the subcarriers unutilized by the 
information transmitter, we can further enlarge the achievable 
rate-energy region. 

We note here that, the result of a one-subcarrier transmission 
strategy of an OFDM system can be regarded as the rank one 
beamforming of the MIMO system in 0. However, the rank 
one beamforming of the MIMO system is in a continuous 
space domain while the subcarrier transmission strategy in the 
OFDM system is in a discrete frequency domain. This results 
in the selection of the optimal subcarrier and the computation 
of the achievable rate-energy (R-E) region being simpler than 
the design of the Geodesic beamformer and computation of 
the R-E region in MIMO 0. 

Organization: Section II presents the system model and 
section III discusses the transmission strategy. Section IV 
presents the simulation results and section V concludes the 
paper. 

Notations: Bold capital letter and lower case letter repre¬ 
sent matrix and vector, respectively. (A)^, tr(A), det(A) 
and IIAII represent the conjugate transpose, the trace, the 
determinant and the 2-norm of a matrix A, respectively. |a| 
denotes the absolute value of a , ( 0 )+ = max(a, 0) and A ^ 0 
indicates that matrix A is a positive semi-definite matrix. The 
matrix with diag(ai,.., oat) is a block diagonal matrix with 
elements a„ and I^r is A x A identity matrix. 

II. System Model 

In this section, we introduce the OFDM system model 
for the SISO two-user IFC where information and energy 
transmitters transmit their OFDM signal to their corresponding 
(information decoding or energy harvesting) receivers. Each 
OFDM symbol has N subcarriers and each transmitter and 
receiver has a single antenna. We assume that each transmitter 
has perfect knowledge of the local CSI, i.e. CSI of the links 
between itself and all receivers. The received signal at i-th 
receiver, G can be written as 

2 

fc=i 

where G is a complex white Gaussian noise vector 

with a covariance matrix For simplicity, we assume 

that (T^ = 1. In addition, we assume that the channel is slow 
fading, the transmitters are perfectly frequency synchronous 
and the cyclic prefix is long enough, such that the channel 


encompassing all subcarriers can be represented by a diagonal 
matrix structure as Hife=diag(/iife^i,Here, hik^n is 
the channel coefficient for the n-th subcarrier between the k- 
th transmitter and the i-th receiver. The relative path loss of 
the cross link compared to the direct link is defined as 5ik G 
[0,1]. The vector G is the transmit signal at the 

k-th transmitter with a transmit power constraint P, such that 


E 


ixfeir 


< P for k = 1, 2. In addition, Xk n is the transmit 


signal for the n-th subcarrier at the k-th transmitter with a 
transmit power pk,n = E . 

When the receiver operates in the ID mode, the achievable 
rate at the i-th ID receiver is given by 


= log2 det(l7v + 


H\ 


( 2 ) 


where = diag(pi i, ..jPi^N) is the power allocation at the 
i-th transmitter with tr(Pi) < P and R_i is the covariance 
matrix of the noise plus interference at the i-th receiver given 

by R_, =lAr + ELiH,fePfcH^. 

i^k 

When the receiver operates in the EH mode, the total 
harvested energy at the i-th receiver (specifically, the total 
harvested energy of the baseband signal as in [l]-[8]) is given 
by 

2 

= C.^tr(H,fcPfcH^-fI^). (3) 

fc=i 

For simplicity, we assume that the energy harvesting efficiency 
constant Ci = 1 the noise power is negligible compared 
to the energy harvested from the EH transmitter such that 

2 

E;,f«^tr(H,fcPfeH^). (4) 

In this paper, we assume that the interference signal from 
other transmitters is not decodable Q. Accordingly, it degrades 
the achievable rate performance at ID user. In contrast, EH 
circuit can harvest energy from the interference signal. There¬ 
fore, the interference signal is beneficial to increase the total 
harvested energy at EH receiver. 

HI. Transmission Strategy 

In this section, we discuss the transmission strategy. There 
are three different scenarios, namely two-ID users, two-EH 
users and one-ID user and one-EH user. In this paper, our 
main focus is on one-ID user and one-EH user scenario. 

A. Necessary Condition For The Optimal Transmission Strat¬ 
egy 

In the two-ID user scenario, the objective is to obtain the 
maximum achievable sum rate subject to the transmit power 
constraint at each transmitter. Since the two users decode data 
information, the total harvested energy is equal to zero. The 
problem in the two-ID user interference channel has been 
considered in several researches. For example, non-cooperative 
game theory has been developed to maximize the achievable 
sum rate with no CSI sharing among transmitters. Each user 


E^ = C.E 


|yiir 
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Fig. 1. Two-user OFDM IFC in {EHi^ID 2 ) mode. 


iteratively updates their power allocation using waterfilling 
algorithm until they reach a Nash Equilibrium ||9|. 

In the two-EH user scenario, the objective is to maximize 
the total harvested energy subject to the transmit power 
constraint at each transmitter. Since the two users harvest 
RE energy, the achievable rate is equal to zero. The optimal 
transmission strategy in two-EH user IFC is for the k-th EH 
transmitter to allocate its full transmit power to the subcarrier 
corresponding to the arg max {|/iife,nP + |/i 2 fc,nP}- 

Now, let us consider one-EH user and one-ID user scenario 
(with their corresponding transmitter respectively denoted as 
EH transmitter and ID transmitter) as in Fig. 1. Without loss 
of generality, we assume that the first receiver is in EH mode 
and the second receiver is in ID mode, (EHi, ID 2 ). The 
discussion is straightforwardly extended to (EH 2 ,IDi). The 
achievable rate-energy region is then given by 

C^R-E){P) ={R< log 2 det (Ijv +H 22 P 2 (R- 2 )-'Hf 2 ) , 

2 

tr (HlfcPfeHf],) > E, tr (P^) < P, 0,k = 1, 2}. 

k=l 

(5) 


Power allocation Pi is coupled with power allocation P 2 
which makes the problem above non-convex with respect 
to Pi and P 2 . The following proposition gives a necessary 
condition for the optimal transmission strategy for the EH 
transmitter at high SNR^. 

Proposition 1: The optimal transmission strategy at high 

SNR for the EH transmitter is to transmit its signal by 

allocating jts transmit power,on,a single subcama'. „ , 
Proof: Let us consider the boundary point ((it, E) of the 

achievable rate-energy, where tr(Hi 2 P 2 Hf^) > E. Then, for 

any given P 2 on the boundary point. Pi = 0 is the optimal 

power allocation at the EH transmitter. That is, because the 

harvested energy from the ID transmitter is sufficient to meet 

the energy constraint E, EH transmitter does not transmit any 

signal causing the interference to the ID receiver. 

For ti-(Hi 2 P 2 Hf 2 ) < E, by letting HaiPiHfi = 
diag(ai,.., on) and H 22 P 2 H|^ = diag(&i,.., 6jv), the bound¬ 
ary point of the achievable rate can be written as 


N 


R= iog2 n (^+ 


2=1 


1 di 


(6) 


'if it is decodable, the interference term can simply be set to 0. 


^Note that in practice, wireless power transfer operates at high SNR [1] 














with the harvested energy from the EH transmitter is given 
by = tr(HnPiHfi) = E- tr(Hi 2 P 2 Hf 2 ). Here, h 
and ai is the received power from ID transmitter and the 
interference power from EH transmitter, respectively. Without 
loss of generality, let us first consider that the EH transmitter 
allocates its transmit power on m subcarriers such that m < N 
and the ID transmitter allocates its transmit power on N 
subcarriers. Eor notation brevity, we sort the subcarrier indices 
at both transmitters such that the EH transmitter allocates its 
transmit power to the hrst m subcarriers. Then, the achievable 
rate in (6) can be written as 

m , , N 

n (i+o.))' <’) 

2=1 * j—m +1 

Because 6^ oc P for i = 1, ■.,N and for a large P, log(l + 
bi) ~ bi, (7) can be written as 

m , N 

2=1 * j—m +1 


In addition, because bi rx P and with a large P, at is negligible 
with respect to bi when E is hnite, then (8) can be written as 


R 



Uhb^ \ 

nr=i(i+«o^' 


(9) 


That is, 01 ^ 1(1 + ill '^^e denominator part of (9) has a 
minimum value when m = 1 and (i.e., oi) is small for 
given E. ■ 

Remark 1; Note that when E is large, both bi and Ui are 
proportional to P. Eor large P, « 1. Then, (8) can be 
written as 

m j N 

p^iog2(n(i+Y^) n (10) 

^ • 1 ^ \ Uji . / 

1—1 j—m-\-l 

where 

/P^\ 

R°^'^og2i^—j={N-m)\og2P. ( 11 ) 


That is, at high SNR and large E, the achievable rate scales 
proportionally with (N — m). Then, it is also maximized when 
TO = 1. 


B. Achievable Rate-Energy Region For One-ID User And One- 
EH User Scenario 


In this section, following the necessary condition of a single 
subcarrier transmission strategy at the EH transmitter, we first 
discuss a single subcarrier selection strategies for the EH 
transmitter and accordingly we propose an iterative algorithm 
to identify the achievable rate-energy region. 

Erom an energy harvesting perspective, maximum harvested 
energy can be achieved when the EH transmitter allocates its 
transmit power to the subcarrier corresponding to the maxi¬ 
mum channel gain (MaxCG) from the EH transmitter to the 
EH receiver. However, this may introduce a large interference 
to the ID receiver. Erom an information decoding perspective, 
the EH transmitter should allocate its transmit power to 


the subcarrier corresponding to the minimum channel gain 
(MinCG) from the EH transmitter to the ID receiver. However 
with MinCG, the EH receiver may not harvest sufficient energy 
to satisfy the energy constraint. In addition, motivated from m 
and 0, the subcarrier index can be efficiently chosen such that 
the Signal to Leakage and Energy Ratio (SEER) is maximized, 
where SEER of the nth subcarrier is defined as 


SLER„= 


\h 


11,22! 




( 12 ) 


Here, ,8 is a fixed regularization parameter. Therefore, the 
subcarrier index n is then chosen as 


(13) 


With the knowledge of the local CSI at the transmitter, we 
then have the following proposition. 

Proposition 2: The optimal fi for a single subcarrier trans¬ 
mission strategy at the EH transmitter is given by 


arg max /in „ ^ 

22 = 1 ,..AT 

for MaxCG 

arg min \h 2 i,n\'^ 

n—l,.,N 

for MinCG 

arg max SLER^ 

n=l,.,Ar 

for SEER 


_ |(^ll,n 

n = arg max- 

|(l21,r2| 


2 ’ 


(14) 


where N is a set of subcarrier indices at the EH transmitter 
satisfying the energy constraint in (5). Here, if N = {}, then 

of the Geodesic energy beamforming in (|[3l. Proposition 3). 
Due to the limited space and the analogy with the proof of 

a new onfl 

subcarrier selection method, namely Signal to Leakage Ratio 
with Energy Constraint (SLREC). n for SLREC is chosen as; 


1) According to the energy constraint in (5), possible 
subcarriers at the EH transmitter are chosen such that 
P\hii^n'^ > E—E 12 , where E 12 is the harvested energy 
from the ID transmitter. 

2) Among the possible subcarriers, n is the subcarrier index 


with the maximum ratio of 


1 ^ 11 ,2t| 

jh^ 


Note that in the SLREC, the searching process is required 
such that a set of subcarriers satisfying the energy constraint 
E is conhned. In contrast, in the SEER, the searching process 
is not required. However, the subcarrier maximizing SEER is 
still implicitly considering E because the required energy E 
is included in the objective (12) of SEER. 

Accordingly, we propose an iterative algorithm that op¬ 
timizes the power allocation at the EH transmitter and the 
ID transmitter iteratively. We have the following optimization 
problem for the achievable rate-energy region of (5) 


(PI) maximize log 2 det {In + H22P2Rl2Hf2) 

2 

subject to ^tr(HifcPfeHffc) >P (15) 

k^l 

tr(Pfc) < -P, PkhO for k = 1, 2, 










where E is the energy constraint and R _2 = lAr+H 2 ir'i±± 2 i- 


iPiHf 

Because the objective function in optimization problem (PI) 
is monotonically decreasing with respect to Pi, we optimize 
Pi using the steepest descent method. If the total harvested 
energy is larger than E, EH transmitter reduces its transmit 
power to minimize the interference to the ID receiver. Given a 


^ — 1 dA.{x) 
dx 


. Then, the 


matrix A(x), ^ IndetA (x) = tr ^A(x) 
gradient of the achievable rate in (15) with respect to is 
given by 

1 / 2 2 — ^ 
\7(Pl,n; P2,n ) — 1 TpT ( Pl,n|^21,n| “f P2,n|^22,n| ) 


ln( 2 ) 

“ (1 + Pl,n|^ 21 ,; 


1^2 


(16) 

Given Pi, the optimization problem (PI) is convex with 
respect to P 2 . Therefore, P 2 at the i-th iteration can be 
solved using the Lagrangian duality method. Given Pi, the 
Lagrangian of optimization problem (PI) is given by 

£(P 2 , A, p) = log 2 det {In + H 22 P 2 (R- 2 )-'Hf 2 ' 
p (tr(P 2 ) - P) + A (tr (Hi 2 P 2 Hf 2 - {E - Pn))) , 
where p and A are the Lagrangian multipliers for the transmit 
power constraint and the energy constraint, respectively. P 2 


(17) 


can be obtained by setting 


dCtPa.A./i) _ 

dP 2 


0 , such that 

+ 


P2= (7-R-2(H22Hf2) ') , 


(18) 


where 7 = (ln(2)(/i/Ar — AHi 2 Hf^)) . Here, /r and A can 
be solved using the subgradient-based method ID, where the 
subgradient of p(A,p) is given by (tr (Hi 2 P 2 H(^) — {E — 
Ell), P — tr(P 2 ))- Note that if the harvested energy from the 
ID transmitter is sufficient, A = 0. Then, P 2 is solved using a 
single-user waterfilling algorithm, maximizing the achievable 
rate with transmit power constraint. 

The iterative algorithm to identify the achievable rate-energy 
(R-E) region is given in Algorithm 1. 

Algo. 1. Identification of the R-E region 

1) Initialize pj°^ such that = P and = 0. 


2) Eor z = 0 : Imax, 
of iterations. 


where is the maximum number 


a) Eor given R^ 2 ’ update P 2 ^ as in (18). 

b) IfELiti-(HiicPi*^Hf,)>P 

Update yJn as in (16). If \yJn = 0, then set 
SjJn = a, where a is a fixed negative parameter. 
Update such that 


3) Einally, 
termined 


H22P^'^(R^^2)“'Hf2),ELitr(HifcPi*^H^ 


= min(max -f A v o) , P). 
achievable 

{R2,Ei) = (iog2ueni7v 

H\ v^2 


the 


as 


rate-energy is 
= (log 2 det (I 

,(i) 

The step size A is given by a fixed value in [0, A^ 
maximum allowable step size, Amax is given by || 2 l 


de- 


E - tr(HnPf^Hfi) - tr(Hi 2 P^^^Hf,) 


(dtj/r' 


AjY,n^X. - 


. The 


(19) 


Due to the discrete subcarrier transmission strategy in 
OEDM system, when EH transmitter transmits its signal at 
the unutilized subcarrier of the ID user, the interference caused 
by the EH transmitter does not create interference at the ID 
receiver. Therefore in step 2. b) of Algorithm 1, when the ID 
transmitter does not allocate its transmit power at subcarrier n 
due to its channel condition, then Jn = 0. Accordingly, the 
EH transmitter will decrease its transmit power at subcarrier 
n by setting yJn = a. 

Algorithm 1 always has negative gradient and positive step 
size, such that Pi is monotonically decreasing. The objective 
function in the optimization problem (PI) is concave with 
respect to P 2 and monotonically decreasing with respect to 
Pi. That is, the objective function of (PI) is quasiconcave 
and the constraints are convex. Therefore, we can conclude 
that Algorithm 1 always converges to a global optimal solution 
for given local CSIT and subcarrier selection strategy ||T]. 

In the simulation results in section IV, an upper-bound on 
the R-E region achievable with the non-cooperative strategy 
with local CSIT and a single subcarrier selection strategy is 
obtained by conducting an Exhaustive Search over all possible 
subcarriers. Specifically, the Exhaustive Search algorithm per¬ 
forms a full search for the optimal single-subcarrier selection 
and power allocation for the EH transmitter, in which the ID 
transmission strategy (i.e. ID transmitter’s power allocation) is 
known at the energy transmitter. Eor a given energy constraint, 
we evaluate the achievable rate from Algorithm 1 for each 
subcarriers at the EH transmitter. Einally, the EH transmitter 
allocates its transmit power at the subcarrier with the high¬ 
est achievable rate. Note that Exhaustive Search algorithm 
requires transmitter cooperation with global CSIT knowledge 
(CSI of the links between all transmitters and all receivers). 

Remark 2: In a single-user waterfilling algorithm, the ID 
transmitter does not allocate its transmit power to the sub¬ 
carrier with poor channel gain-to-noise ratio depending on 
the water level. That is, if there is some interference on 
the subcarriers unutilized by the ID transceiver, it will not 
degrade the achievable rate performance. Hence, \/Jn = 0. 
However, for the EH transmitter to make use of the knowledge 
of those subcarriers unutilized by the ID transmitter, the 
EH transmitter requires some transmitter cooperation (as in 
the Exhausive Search algorithm) where the indices of those 
unutilized subcarriers are informed to the EH transmitter. 

C. Information Sharing of Unutilized Subcarrier Indices 

Erom Remark 2, by allowing the ID transmitter to share with 
the EH transmitter the indices of the subcarriers unutilized by 
the ID transmitter, we can further enlarge the achievable R-E 
region. Let Z be the set of unutilized subcarrier indices at the 
ID-user, i.e., Z = {n„|P 2 .n„ = 0}. R in (5) can be written as 

N 

R = X! (1 +P2.n|^22,nP(l +Pl,n|/l21,nn (20) 

n^'Zi 

That is, the achievable rate is not degraded by the interference 
signal on the subcarrier unutilized by the ID transceiver. There¬ 
fore, if the EH transmitter knows the unutilized subcarrier 







R-E tradeoff curves for one-ID user and one-EH user 


indices (i.e., Z), it can transfer the energy without interfering 
the ID operation. Accordingly, with a transfer of information 
from the ID transmitter to the EH transmitter, we can allocate 
the power to the subcarrier unutilized at the ID transceiver. 
The subcarrier index n for SLREC with information sharing is 
chosen as: According to the energy constraint in (5), possible 
subcarriers at the EH transmitter are chosen. Among them, 

1) If there are unutilized subcarriers at the ID transmitter, 
choose n as the subcarrier index with the largest channel 
gain 

2) Else, choose n among them as the subcarrier index with 

the maximum ratio of (i. e., SLREC). 

|fl21,n| 

The iterative algorithm to identify the achievable rate-energy 
region for SLREC with information sharing of unutilized 
subcarrier indices at ID transmitter is given in Algorithm 1. 
IV. Simulation Results 

In this section, we present the simulation results. We assume 
that the path loss is set to such that the channel = 

where Hife = diag(F[hife On-l]"^)- hik G 
is the multipath channel of length L = 3 where the 
frequency selective fading channel is modeled using three- 
tap exponentially distributed power profile, each with complex 
zero-mean random Gaussian distribution. F S 
normalized Fourier matrix. Here, 5ik — 1 for i = k and 5ik 
— 0.8 for * ^ fc. In addition, P is set to SOmVE, the noise 
power at each subcarrier is set to IpW and N is set to 8. 

Fig. 2 shows the R-E tradeoff in the two-user OFDM IFC. 
The EH transmitter allocates its transmit power on a single 
subcarrier based on MaxCG, MinCG, SEER and SLREC 
as described in section III-B and the achievable rate-energy 
region is determined based on Algorithm 1. The simulation 
result shows that the MinCG increases the achievable rate at 
the ID receiver while the MaxCG achieves a larger harvested 
energy at the EH receiver. As expected, the R-E region for 
the SEER covers both MaxCG and MinCG R-E regions. The 
proposed SLREC subcarrier selection method has a larger R-E 
region than SEER. That is, SLREC minimizes the interference 
gain at ID receiver provided that the energy constraint is 
satisfied. In addition, in the region where the harvested energy 
is less than 18/iW, the achievable rate is unchanged. Here, be¬ 
cause the harvested energy from the ID transmitter is sufficient 
to meet the energy constraint, EH transmitter does not transmit 
any signal to prevent any interference to ID receiver. 

Fig. 3 shows that SLREC with information sharing of the 
unutilized subcarriers at the ID transmitter has a larger achiev¬ 
able R-E region than SLREC without information sharing. 
That is, because the interference signal at the unutilized sub¬ 
carriers at ID-user does not degrade the achievable rate, the EH 
transmitter can transfer the energy without compromising the 
achievable rate at the ID receiver. In addition, the Exhaustive 
Search algorithm provides an upper-bound on the achievable 
R-E region with a single subcarrier selection strategy at the 
EH transmitter. 


V. Conclusion 


In this paper, we have discussed SWIFT in SISO two- 
user OFDM non-cooperative interference channel with perfect 



Fig. 2. R-E tradeoff for MaxCG, MinCG, SLER and SLREC. 


R-E tradeoff curves for one-ID user and one-EH user 



Fig. 3. R-E tradeoff for SLREC, Exhaustive Search and SLREC with 
information sharing by ID-user. 


knowledge at the transmitters of the (local) CSI of the link 
between itself and all receivers. We have found that, the 
necessary condition for the optimal transmission strategy in 
high SNR is for the EH transmitter to allocate its transmit 
power on a single subcarrier. Accordingly, we identify the 
achievable rate-energy regions for different subcarrier selec¬ 
tion strategies - MaxCG, MinCG, SLER and SLREC. We have 
found that SLREC exhibits higher R-E performance than the 
other subcarrier selection strategies. Interestingly, with some 
transmitter cooperation allowing the indices of the subcarriers 
unutilized by the ID transmitter to be shared with the EH 
transmitter, we can enlarge the achievable rate-energy region. 
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